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ABSTRACT 

We present VLT-UVES echelle spectroscopy of the Hi and C IV absorption in the spatially- 
extended Lya emission around two high-redshift radio galaxies 0200H-015 (z=2.23) and 0943- 
242 (z=2.92). 

The absorbers in 0943-242 exhibit little additional structure compared with previous low- 
resolution spectroscopy and the main absorber is still consistent with Hi column density of 
^ 10^^ cm^^ . This is consistent with a picture in which the absorbing gas has low density and 
low metallicity and is distributed in a smooth absorbing shell located beyond the emission- 
line gas. However, the main absorbers in 0200H-015 are very different. The previous single 
absorber fit of Hi column density ~ 10^^ cm^^ , now splits into two ~ 4 x 10^'' cm~^ ab- 
sorbers which extend more than 15 kpc to obscure additional Lya emission coincident with a 
radio lobe in these high-resolution observations. Although consistent with the shell-like dis- 
tribution for the absorption systems, 0200H-015 requires a much higher metal enrichment than 
0943-242. The metallicity, inferred from the C IV absorption, is considerably lower in 0943- 
242 than in 0200H-015. We explain these differences with an evolutionary scenario based on 
the size of the radio source. In both sources the Hi absorption gas originates from either a 
gas-rich merger or pristine cluster gas which cools and collapses towards the centre of the 
dark matter halo. The higher metallicity in the larger radio source (0200H-015) may be a result 
of a starburst driven superwind (concurrent with the triggering of the radio emission) which 
has engulfed the outer halo in this older source. 

We also find a significant blue asymmetry in the HeIIA164G emission line, suggesting 
that the Une emitting gas is outflowing from the central regions. Dust obscuration toward the 
central engine, presumably due to the dusty torus invoked in Unified Scheme, prevents us 
from seeing outflow away from our line-of-sight. 

Key words: galaxies: active - galaxies: halos - galaxies: high-redshift - galaxies: emission 
lines - galaxies - absorption lines 



1 INTRODUCTION 

The existence of powerful radio galaxies at high redshift (z > 2) 
demonstrates that a population of supermassive black holes (> 
10^ Mq ; Dunlop et al. 2002) was in place just a few billion years 
after the Big Bang. Like their counterparts at low redshift, the high- 
redshift radio galaxies (HzRGs) appear to reside in the most mas- 
sive elliptical galaxies at their epoch (e.g. Jarvis et al. 2001b), 
in line with the established correlation between black hole and 
spheroid mass (Magorrian et al. 1998). Whilst the discovery of ra- 
dio quiet galaxies at high redshift via Lyman drop-out techniques 
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means that they are no longer our only probe of high redshift galaxy 
formation, the HzRGs still provide the most important insights into 
the early formation of the most massive bound structures. Indeed, 
the depth of the gravitational potential wells in which they form 
renders them prime targets for searching for high-redshift proto- 
clusters, as confirmed by the recent discoveries of over-densities 
of LyQ emitters around such objects (see e.g. Kurk et al. 2000; 
Pentericci et al. 2000; Venemans et al. 2002). 

One of the most prominent characteristics of HzRGs are 
their extended emission line regions (EELRs), which are luminous 
(> 10'^'^ W in Lyof), tens to several hundred kpc in size (often 
aligned with the radio axis) and kinematically active (FWHM > 
1000 km s~^). Their dominant ionisation mechanism seems to 
change from shocks to photoionisation by the central engine as 
the radio source expands beyond the confines of the host galaxy 
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(Best et al. 2000; De Breuck et al. 2000; Jarvis et al. 2001a). Like- 
wise the kinematics of the gas may reflect the competing influences 
of gravity, and energy input via shocks from the radio source and 
starburst-driven superwinds. However, there are many unresolved 
issues which remain concerning the structure, origin and fate of the 
emission line gas, which may hold important clues to many aspects 
of massive galaxy formation. Did the gas originate in a cooling flow 
or a merger, or was it expelled from the central galaxy during a vio- 
lent starburst? Is it in the form of cloudlets, filaments or expanding 
shells of material, and what is their composition? What is the ulti- 
mate fate of the gas, will it form stars or escape from the galaxy to 
enrich the surrounding intracluster and intergalactic media? 

A new perspective on many of these issues was opened up with 
the discovery by Rottgering et al. (1995) and van Ojik et al. (1997; 
hereafter v097) that most of the smaller high redshift radio sources 
(those with projected linear size D < 50 kpc) exhibit spatially re- 
solved HI absorption in their Lya emission-line profiles, with col- 
umn densities in the range lO^^-lO^'' '^ cm"^ . Binette et al. (2000; 
hereafter BOO) found, for the z = 2.92 radio galaxy 0943-242, 
C IVAA1548,1551 absorption lines superimposed on the C IV emis- 
sion, at the same redshift as the main Lya absorption system. 
They could not reconcile the observed Civ/Lya emission- and 
absorption-line ratios with a model in which the absorption- and 
emission-line gas are co-spatial. Instead they proposed that the ab- 
sorbing gas is of lower metallicity (Z ~ O.OIZq) and located 
further away from the host galaxy than the emission line gas, be- 
yond the high pressure radio source cocoon. This material, they 
claimed, is thus a relic reservoir of low metallicity, low density 
(~ 10~^'^ cm~^ ) gas, similar to that from which the parent galaxy 
may have formed. This shell-like distribution of the absorbing gas 
was also invoked by Jarvis et al. (2001a) and De Breuck et al. 
(2000) to explain various correlations found in a complete sample 
of high-redshift (z > 2) radio galaxies. 

In this paper we probe the HI and C IV absorbers in greater 
detail using the UVES spectrograph on the VLT to obtain echelle 
spectra of two high-redshift radio galaxies, the aforementioned 
0943-242 and 0200-1-015 at z=2.23. The new observations offer a 
factor of 10 improvement in spectral resolution over previous in- 
vestigations and thus have the potential (i) to reveal whether the 
main absorbers are genuinely in the form of a single component 
with ordered global kinematics, or whether they fragment into a 
number of weaker absorbers; (ii) to probe the emission-line kine- 
matics with high velocity resolution. 

In section ^ we describe the observations and data reduction 
of our VLT-UVES spectroscopic observations and in section ^ we 
briefly discuss the emission- and absorption-line fitting procedure. 
In section ^ we use our new high-resolution spectra to investigate 
the absorption halo of the z = 2.92 radio galaxy 0943-242. Sec- 
tion ^ is a detailed discussion of the absorbing halo around the 
z = 2.23 radio galaxy 0200-1-015 and we also investigate the kine- 
matics of the narrow-emission lines via the unabsorbed Hell emis- 
sion line and briefly discuss the effect of emission-line asymmetry 
on our absorption line fitting. The origin and fate of the absorbing 
halos is discussed in section ^ and we provide a summary of our 
conclusions in section ^. 

All physical distances are calculated assuming Ho = 
70kms"^ Mpc"\nM = 0.3 and J^a =0.7. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were performed with the UVES echelle spec- 
trograph on VLT UT2 at the European Southern Observatory on 
the night of 2001 December 8-9. For 0200-1-015, dichroic 1 was 
used with cross-disperser 1 in the blue arm (central wavelength 
of 3950A, for Lya), and cross-disperser 3 in the red arm (central 
wavelength of 5800A, for C IVAAI548, 155 1 and HellA 1640); three 
separate exposures of I hour were taken. For the 0943-242, just the 
red arm was used, with cross-disperser 3 set to a central wavelength 
of 5200A, sufficient to cover both Lya and C IVAA 1548, 1551. Sep- 
arate exposures totalling 3.4 hours were acquired. For all observa- 
tions on-chip binning of 2 x 3 (spatial x spectral) was used, re- 
sulting in a pixel size of 0.5 arcsec in the red arm, and 0.36 arc- 
sec in the blue. The resulting pixel scales in the dispersion direc- 
tion for the various set-ups were in the range 0.05-0.06A. Together 
with a slit- width of 1.2 arcsec, the resulting spectral resolutions 
were 25000-40000. The seeing was steady throughout the night at 
~ 0.8 arcsec, and the sky dark and photometric. The slit position 
angles were chosen to align with the radio axes and previous obser- 
vations, namely 155 degrees for 0200-1-015 and 74 degrees east of 
north for 0200-1-015 and 0943-242 respectively. 

Offline reduction was performed with IRAF, following the 
procedures for echelle data described by Churchill (1995). For each 
target, the individual exposures were median combined to remove 
cosmic rays, bias-subtracted and flat-fielded. The order-definition 
frames were used to determine the locations of the orders, which 
were then extracted and wavelength calibrated using a ThAr arc 
and combined. The blaze function of the grating was removed using 
observations of the G-type spectrophotometric standard LTT 1020. 
Spectra for individual spatial regions along the slit were extracted, 
rebinned by a factor of 2 in dispersion to improve the signal-to- 
noise ratio. 



3 ANALYSIS 

The extracted 1-d spectra were fitted with a series of Voigt absorp- 
tion line profiles superimposed upon a Gaussian emission envelope. 
We choose an underlying Gaussian emission-line profile as studies 
of narrow-line profiles in Seyfert-2 galaxies show that this is pre- 
ferred over a Lorentzian profile in the majority of cases (e.g. Veron- 
Cetty, Veron & Gon^alves 2001), but see section for further dis- 
cussion of this point. The model profiles were convolved with the 
instrumental function prior to fitting (as described by Rottgering et 
al. 1995 and BOO). There is no detected continuum underlying the 
emission lines. 



4 0943-242: CONCLUSIVE EVIDENCE FOR A DISTINCT 
SHELL OF ABSORPTION-LINE GAS? 

Fig. [l] shows the 1-d spectrum of the Lya and CiV emission 
lines in 0943-242, extracted over the central 5 arcsec of the slit 
which covers all of the emission. The fit, with a series of Voigt pro- 
file absorbers superimposed upon a Gaussian emission envelope is 
shown, and the parameters tabulated in Table nl. The errors were 
calculated by assuming a distribution of Ax (= ~ Xmin) 
(Lampton, Margon & Bowyer 1976). The error on each parameter 
was calculated by setting Ax^ ~ 1 and allowing the other param- 
eters to float. The errors quoted are la. 

The most notable feature of this spectrum is that the main ab- 
sorber remains as a single system of column density ~ lO^'' cm~^ , 
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0943-242 . 



Absorber 


z 


b 


log Af (HI) 






(kms~-^ ) 


(cm~^ ) 


1 


2.9066 ± 0.0062 


88 ±45 


14.02 ± 0.30 


2 


2.9185 ±0.0001 


58 ±3 


19.08 ± 0.06 


3 


2.9261 ± 0.0005 


109 ± 35 


13.55 ±0.16 


4 


2.9324 ± 0.0001 


23 ± 17 


13.35 ±0.30 








log Af (CIV) 








(cm~^ ) 


2 


2.9192 ±0.0001 


119±2 


14.58 ± ± 0.04 



Table 1. Parameters of the Voigt absorption profile fits for the main com- 
ponent of 0943-242. The final row con'esponds to the fit of one absorber to 
the CIV profile. 

and is completely black at its base with no evidence for substruc- 
ture. No new absorption systems are identified in addition to those 
found by Rottgering et al. (1995) at a factor of ~ 10 lower resolu- 
tion. 

This provides strong evidence that these absorption systems 
are physically distinct from the extended narrow-emission-line re- 
gion. If the absorption was caused by gas which is mixed with the 
emission-line clouds, as postulated by v097 then we would expect 
a series of narrow-absorption troughs at various wavelengths across 
the Lya emission profile. The fact that we see one main absorber, 
with no detectable substructure, shifted blueward of the emission- 
line peak provides compelling evidence that the absorption gas en- 
compasses the whole of the emission-line regions with a covering 
factor of unity along our line-of-sight. This is in line with the argu- 
ment proposed by BOO who fitted the emission- and absorption-line 
profile of the C IVAA1548,1551 doublet in 0943-242 to determine 
the metallicity of the absorption systems compared to the emission- 
line gas. Their results show that the emission and absorption line 
ratios of C IV and Lya are incompatible with photoionisation or 
coUisional ionisation of cloudlets with uniform properties and the 
possibility that the absorption and emission phases are co-spatial is 
rejected. A different model was preferred in which the absorption 
gas has low metallicity and is located further away from the host 
galaxy than the emission line gas. 



5 0200+015: A DETAILED STUDY 

In this section we provide an in depth study of the absorption sys- 
tems around the z — 2.23 radio galaxy 0200-1-015. This object was 
taken from the original sample of v097 and was chosen on ac- 
count of its interesting absorption-line profile in the ^ 2A reso- 
lution spectrum of v097. A further advantage of using this object 
for our study is that its redshift permits further observations of all 
of the principal optical emission lines in the near-infrared transmis- 
sion windows arising from AGN phenomena which will enable us 
to determine line ratios and investigate the ionisation mechanisms 
at work as a function of radial distance from the central source. 

Examination of the Lya region of the CCD frame of this ob- 
ject shows that the emission deviates from the centre of the order 
(as defined by the order definition frame) as a function of wave- 
length. This reflects spatial extent in the emitting and absorbing 
gas. Spectra were extracted for 2 spatial regions: the main body of 
the galaxy between ±1.25 arcsec of the centroid, and a secondary 
component lying between 1.25-2.5 arcsec to the north-west (1 arc- 
sec translates into 8kpc at this redshift), coincident with one of 
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Figure 1. (top) The Lya profile of 0943-242, with the absorption 
model overlaid, (bottom) The CIV profile of 0943-242, with the 
absorption model overlaid. 



Table 2. Parameters of the Voigt absorption profile fits for the main compo- 
nent of 0200-hO 15 



Absorber 


z 


b 


log V(HI) 






(kms^^ ) 


(cm^^ ) 


1 


2.2239 ± 0.0007 


612 ±70 


14.96 ± 0.06 


2 


2.2282 ± 0.0002 


116 ±27 


14.73 ±0.14 


3 


2.2307 ± 0.0001 


61 ±8 


1458 ±0.16 


4 


2.2349 ± 0.0001 


86± 11 


14.38 ±0.10 



the two radio lobes (see the radio continuum image in Carilli et 
al. 1997). Spectra with overlaid absorption models for these two 
components are shown in Fig. ^ and the parameters tabulated in 
Tables I and I 

5.1 The Lya emission from the nucleus 

The core Lya emission-line profile (i.e. the emission asso- 
ciated with the centre of the radio galaxy host) has a red- 
shift of 2 = 2.2292 and full width half maximum FWHM = 
1490 kms~^ . The emission line exhibits a number of absorbers 
within 320 km s~ ^ of the redshift of the peak of the Lya emission, 
consistent with the lower resolution spectrum of v097. However, 
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Figure 2. (a) (top) The Lya profile of the central 2.5 arcsec of 
0200+015, with the Gaussian emission envelope and the absorption 
model overlaid, (b) (bottom) As in (a) but for the secondary (offset) 
component of 0200+015, with the absorption model overlaid. 



at this higher resolution we do find that the main absorber, with a 
measured column density of A^'hi ~ 10^^'^ cm~^ from v097, is ac- 
tually comprised of two absorption troughs with column densities 
of Nm = 10" ''^ cm~^ and iVni = IQ^-^s cm"^ with velocity 
widths of 116 km and61kms"^ and redshifts of 2: = 2.2282 
and z — 2.2307 respectively. 

There is also a A^hi = 10^* cm~^ absorber blueward of these 
lines if our Gaussian emission profile is accurate. This absorption 
trough has a much higher Doppler parameter of 6 = 612kms~^ , 
although the sensitivity of this absorption profile on the fit to the 
emission line is significant and we explore this further in sec- 
tion O 



5.2 Emission-line asymmetry and the evidence for outflow 

One of the very noticeable aspects of the spectroscopic data on 
these two high-redshift radio galaxies is that in the Hell emission 
line, in which associated absorption does not affect the profile of 
the emission line, the blue side of the emission line has a longer 
tale than the redward side (Fig. Using the profile analysis origi- 
nally defined by Heckman et al. (1981) we assign a relative asym- 
metry, AI, to the emission at varying fractions of the total inten- 
sity where AI = 2[Ac — Xc{i)]/ijj{i). Here Ac(i) and Lu{i) are the 
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Figure 3. The Hell profile of the central 2.5 arcsec of 0200+015. 



Table 3. Asymmetry of the Hell line in 0200+015 



per cent of peak ((') 


Ac /A 




AI 


10 


5296.4 


43.6 


29.4 


20 


5296.9 


36.1 


32.4 


33 


5298.1 


29.9 


31.8 


50 


5298.9 


23.0 


33.9 



line centre and width at relative intensity i, respectively, and Ac is 
the wavelength of the peak of the profile. The results of the profile 
analysis are presented in Table ^. Using the FWHM of the asym- 
metric profile to determine the velocity dispersion of this line we 
findv^ 1920kms"i . 

This blue-tail line asymmetry has been found previously for 
the narrow-line regions in many Seyfert galaxies (e.g. De Robertis 
& Shaw 1990) and is probably associated with an inflow or outflow 
of the gas into or from the central regions. Evidence from narrow- 
line Seyfert galaxies points towards the view that the asymmetry is 
caused by outflowing gas propagating towards us along our line- 
of-sight. The counter flow, which would produce a red-tail asym- 
metry, is obscured by dust residing within the radius of the narrow- 
emission-line region. This is in line with the dusty torus invoked 
in orientation based Unified Schemes. The alternative view, that of 
inflowing gas, requires the dust to be mixed with the emission-line 
clouds. This would mean that the blue tail emanates from the far 
side of the galaxy. In this case the illuminated clouds on the near 
side along our line of sight suffer higher extinction than the illu- 
minated clouds on the far side (e.g. De Robertis & Shaw 1990). 
However, ISO observations of NGC 4151 show that this explana- 
tion is not sufficient to explain the emission-line characteristic in 
the infrared where the emission lines are much less susceptible to 
obscuration. Sturm et al. (1999) propose a central, optically thick 
obscuring screen close to the nucleus with outflowing emission- 
line regions, most plausibly associated with a nuclear outflow or a 
starburst. 

Obviously, this particular characteristic of the emission line 
region needs exploring further with high-resolution spectroscopy 
of emission lines with a range in ionisation parameters to determine 
whether the velocity of the material is dependent on the density of 
the gas and its distance from the central ionizing source. However, 
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Table 4. Parameters of the Voigt absorption profile fits for 0200+015, as- 
suming asymmetric emission-line profile using the asymmetry parameters 
from the Hell emission fine. The top panel is for the four absorbers in the 
Lya emission line from the nucleus , the lower panel represents the single 
absorber in the CIV profile. 
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Figure 4. The Lya profile of the central 2.5 arcsec of 0200-1-015 
fitted with an asymmetric Gaussian with the same asymmetry as 
the Hell emission line . 



for the purposes of this paper where we are predominantly con- 
cerned with the absorbing shells, it highlights the fact that assuming 
a simple Gaussian profile for the narrow emission lines may be too 
simple. We have therefore refitted the Lya and C IV emission lines 
for our two objects to determine whether this amount of asymme- 
try may significantly affect our fitted parameters to the absorption 
troughs in any way. 

We refitted our emission-line profiles with two single tailed 
Gaussians with common peak intensities and a fixed ratio of 
FWHM derived from the asymmetry in the Hell line as a simple 
approximation to the asymmetry profile. Although we note that the 
different ionisation potential for Lya and C IV compared to Hell 
may affect the asymmetric profile, we do not expect this to be a sig- 
nificant difference given the uncertainties in our model fitting. With 
this fixed asymmetry the resulting fits to the absorption troughs in 
both Lya and C IV are not significantly altered, these fitted param- 
eters are listed in Table |3| 

We therefore conclude that although there is evidence for line 
asymmetry in the narrow-emission-lines in HzRGs this does not al- 
ter the fits to the absorption troughs significantly and the original 
results are still valid. This is quite surprising given the uncertain- 
ties arising from the emission-line profile. However, we note that 
it is very difficult to change the actual absorption columns given 
any reasonable Lya emission-line profile. This is particularly true 
where the emission is at a low level, i.e. around absorber #1 in 
Fig. ^ and Fig. ^. Whereas the largest errors occur near the peak of 
the emission where the absorbing columns (absorbers #2 and #3) 
have larger associated errors. Furthermore, the fits to absorbers #2 
and #3 are more tightly constrained in b because the higher degree 
of absorption at lower b means that there is little room to manou- 
vre in the absorption line fitting. The combination of these effects 
results in little variation in the absorptions line column densities. 
However, we note that more extreme Lya emission-line profiles 
will result in larger errors on the absorption-line fits in the wings of 
the emission line. 



5.3 The Lya emission associated with the radio lobe 

We also observe relatively bright Lya emission at the position of 
one of the radio lobes associated with 0200-1-015. This emission lies 



Absorber 


Line 


z 


b 


log N{m) 








(kms-l ) 


(cm~^ ) 


1 


Lyo 


2.2206 ± 0.0012 


874 ± 130 


14.93 ± 0.07 


2 


Lya 


2.2282 ± 0.0003 


120 ± 25 


14.71 ±0.12 


3 


Lya 


2.2307 ± 0.0001 


55 ±7 


14.52 ±0.17 


4 


hya 


2.2349 ± 0.0001 


83 ± 11 


14.39 ± 0.09 










log Ar(CIV) 










(cm-2) 


2 


CiV 


2.2285 ± <0.0001 


59 ±2 


14.64 ± 0.01 



~ 2.5 arcsec (« 20 kpc) away from the nuclear emission. The best- 
fit Gaussian profile to this emission line shows that the emission has 
a slightly higher redshift than that found towards the central radio 
nucleus {Sv ~ 150 kms"^ ) and has a FWHM = 1580 kms"^ . 

This emission is presumably due to the passage of the radio 
emitting plasma disrupting the gas at these radii. Invoking the sce- 
nario set out by Jarvis et al. (2001a) based on the model of Bremer, 
Fabian & Crawford (1997), in which the very extended narrow- 
emission-line regions are compressed and shredded by the radio 
jets, then this line emission may have been associated with the in- 
nermost low-density absorbing shell, which has condensed enough 
to be seen in emission. This would also explain the slight veloc- 
ity shift of this emission line if the radio jet, which has essentially 
boosted this emission, is propagating away from our line-of-sight. 
However, regardless of the process by which the Lya photons are 
emitted, the surrounding absorbers are distributed in the same way 
as we see in the nuclear Lya emission. This is probably the most 
compelling evidence that the absorbing gas is distributed in a shell 
encompassing the whole of the radio galaxy with a covering factor 
of ~ 1, in line with our conclusions for 0943-242. The velocity 
difference of the absorption systems detected in the Lya emission 
associated with the nuclear emission and that at the position of the 
lobe is ~ 30 kms^^ , thus we can be quite confident that the ab- 
sorption systems seen in the lobe emission are associated with the 
absorption systems associated with the nuclear emission. 

We also find more absorbers in this extended emission line in 
addition to the four main absorbers found in the nuclear Lya emis- 
sion line. These additional absorbers are on the red side of the Lya 
emission line and have significantly lower column densities. The 
largest velocity offset from the Lya emission-line redshift for these 
redshifted absorption systems is ~ 600 km s~^. The lower mea- 
sured column densities of the absorption gas, implies that the cover- 
ing factor is lower along the line-of-sight to the emission associated 
with the lobe. 

If we assume that the emission from the lobe, which is red- 
shifted with respect to the Lya emission-line redshift at the nucleus, 
is on the far side and is propagating away from us, then any emis- 
sion from this position would have to pass through both a lower 
column density redshifted HI component associated with the far 
side and the bulk of the gas on the near side which is blueshifted. 
This would explain the higher column density absorbing gas on the 
blue wing of the emission line. 
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Table 5. Parameters of the Voigt absorption profile fits for the secondary 
component of 0200+015 



Absorber 


2 


b 

(kms-i ) 


log N(m) 
(cm-2) 


1 


2.2241 ± 0.0010 


216 ± 110 


14.68 ± 0.26 


2 


2.2276 ± 0.0002 


106 ± 29 


14.52 ± 0.22 


3 


2.2307 ± 0.0001 


61 ± 11 


14.70 ± 0.24 


4 


2.2347 ± 0.0002 


111 ±21 


14.63 ±0.16 


5 


2.2387 ± 0.0003 


44 ±37 


13.69 ± 0.35 


6 


2.2402 ± 0.0008 


24 ±22 


13.23 ± 0.65 


7 


2.2423 ± 0.0018 


87 ±59 


13.95 ±0.58 
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Figure 5. The CIV profile of 0200+015, with the vertical ticks indi- 
cating the location of the absorption doublet. 



5.4 The CIV emission from the nucleus 

In the red arm of the spectrum, emission lines of the 
ClvAA1548A,155lA doublet and Hell A 1640A are present. The 
former is shown in Fig. ^ extracted across the entire galaxy be- 
cause no spatial structure was evidenl|^ and rebinned in wavelength 
by a factor of 10 to improve the signal to noise ratio. From our line 
fitting we find that the centre of the emission line is redshifted with 
respect to the nuclear Lya emission line by ~ 140 kms~^ , how- 
ever this may be an artifact of trying to fit a symmetric line profile 
(cf. section Isl) with a FWHM = 1940 km s"^ 



A clear absorption doublet can be identified, corresponding to 
ClvAA1548A,155lA at z ~ 2.2285, i.e. roughly corresponding 
to absorption system #2 in Lya (see Tables ^and ^. C IV absorp- 
tion due to Lya absorber #3 would appear at ~ 500lA,5010A but 



^ Using the Lya: C IV emission-line flux ratio at the nucleus, our obser- 
vations are not sensitive enough to detect the expected flux in C IV at the 
position of the radio lobe. 



Table 6. Parameters of the Voigt absolution profile fits for the CIV compo- 
nent of 0200+0 1 5 . 



Absorber 



b log V(CIV) 

(kms^-'- ) (cm^^ ) 



no obvious signature is present. Lyo? absorber #4 would appear 
at ~ 5007A,5017A respectively, thus the 1548A trough would 
overlap on the 1551 A trough of absorber #2. However, the lack of 
any apparent absorption at ~ 5017A suggests that the C IV ab- 
sorption column is significantly lower (and is consistent with no 
C IV absorption at all) in this shell as compared to that associated 
with absorber #2, and we assume that all of the C IV absorption can 
be accounted for by associating it with Lya absorber #2. This im- 
plies a metallicity gradient in the environment of this radio source. 
Overzier et al. (2001) recently found strong evidence for a metal- 
licity gradient within the gaseous halo of the z = 2.49 radio galaxy 
MRC2104-242. Overzier et al. (2001) attributed this metallicity 
gradient to a scenario in which the gas is associated with a mas- 
sive cooling flow or originates from a merging event. This scenario 
can also be extended to the absorption gas, and we return to this 
point in section |^ However, one should keep in mind that the mea- 
sured A'^civ is an upper limit, due to possible extra absorption from 
the other absorption shells. 



5.5 The metallicity of the absorbing shell 

In this section we investigate whether the absorption- and emission- 
line gas can exist co-spatially or whether they are distributed differ- 
ently based on ionisation models of gaseous slabs with the column 
densities we observe. We use the method of BOO to investigate the 
metallicity of the absorbing shell #2 which is the only one clearly 
seen in absorption in C IV. We use the observable quantity F as de- 
scribed by BOO to describe the various ratios between the emission- 
and the absorption-line properties of both Lya and C IV. T is de- 
fined as 

Iciv Nhi 



^ Lya 



Nci 



(1) 



2.2285 ± < 0.0001 69 ± 1 



14.69 ± 0.02 



where Iciv / Ihya is the ratio of the emission-line fluxes and 
Nhi / Nciv is the ratio of the measured absorption columns. 
These four quantities carry information on the three ionisation 
species H^ H+ and C+^ 



5.5.1 Co-spatial absorption and emission line gas 

Using the current spectra and the Gaussian fit to the emission lines, 
we obtain Iciv / Ihya ~ 0.12 ± 0.10^. Adopting the column of 
Lya absorber #2 in Table ^ and that of CIV in Table |^ we obtain 
Nhi/Nciv = 1.09±0.04, hence r 0.13±0.12. This is a factor 
of ~ 5 X 10* smaller than the value found by BOO for 0943-242 of 
r = 5400. BOO used the high P value in 0943-242 as evidence that 
the absorbing gas was physically distinct from the emission line re- 
gions and was of much lower density and metallicity. Calculations 
with the MAPPINGS Ic code (Ferniit et al. 1997) indicate that our 
extremely small value of F effectively rules out that the same so- 
lar metallicities characterize both the emission- and the absorption- 
line regions. In fact, reproducing a F of 0.13 requires a metallicity 
of 10 Zq if we insist that the absorption and emission gas com- 
ponents have the same metallicities. Therefore, it is theoretically 
possible to reproduce F with a strong overabundance of metals, 
which may suggest that the emission gas and absorption gas might 
be more directly inter-related in 0200+015 than in 0943-242. On 
the other hand, the fact that the emission has such a large FWHM 



The error quoted is the Icr uncertainty derived from the symmetric and 
asymmetric Gaussian fits to the emission-lines of both Lya and C IV. 
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as compared to the b values of the absorption profiles and the near 
unity covering factor of the absorption gas, would still favour the 
simple picture in which the absorption gas is more quiescent, pos- 
sibly situated in quasi-uniform shells which are either expanding or 
inf ailing. 

5.5.2 An absorbing shell surrounding the emission-line region 

We now consider the absorption gas in isolation and attempt to ex- 
plain the process by which the C IV column can be approximately 
equal to that of H I in the case of absorber #2. At a given metallic- 
ity, there exists a physical upper limit to the Nciv I Nhi column 
ratio which photoionisation models cannot exceed. This is shown 
in Fig. ^ which show photoionisation calculations by an ionizing 
continuum of photon index a = — 1 (F„ oc u") as a function of 
the ionisation parameter^ U . The geometry considered was that of 
an optically thin slab with all models having the same thickness in 
Nhi of 10^^'^'^ cm^^. Inspection of Fig. n shows that the solar 
metallicity curve cannot reach the observed value of 0.9 while a 
slab of 10 Zq can again reproduce the target ratio. If the ionisation 
process, however, was dominated by coUisional ionisation, then one 
could not rule out solar metallicity. This is shown by the dashed- 
line in Fig. ^ which represents a slab which is initially photoionised 
with a low value of U but where the temperature is imposed. The 
temperatures are increasing (starting from the equilibrium value of 
13000 K) in steps of 0.1 dex and collisional ionisation progressively 
takes over. The open squares identify the first eight models, up to 
8.2 X 10^ K. To show that the column ratio actually decreases at 
much higher temperatures, we used filled triangles to denote the 
models in which the temperature starts from 1.03 x 10^ K and in- 
creases in steps of 0. 1 dex. We have verified that, as expected, the 
peak reached in the column ratio is not higher for a slab with higher 
ionisation parameter (the peak is actually lower for U > 0.01). 

We infer from the above calculations that a solar metallicity 
would therefore require an unusual fine tuning in order to give a 
temperature which lies in the range of ~ 0.7-1.4 x 10^ K. This is 
not a very probable state for any astrophysical plasma since this is 
exactly the same temperature range across which the cooling curve 
is known to peak (e.g. Fig. 8 in Ferruit et al. 1997). For the gas to 
be maintained in this temperature range would require an additional 
(unknown) heating mechanism. Cooling shocks can in principle re- 
produce the observed C IV columns but the models so far computed 
indicate huge Nhi columns ~ 10^® cm~'^ (see Table 2D in Dopita 
& Sutherland 1996), a value which is many orders of magnitude 
beyond the observed value for absorber #2 of lO^*'^'^ cm~^. For 
these reasons, it appears to us more plausible to invoke metal-rich 
absorbing gas. 

Therefore, although the value of T in 0200+015 (unlike 0943- 
242 cf. BOO) can be reconciled with a co-spatial distribution of 
absorption- and emission-line gas, this is not our preferred scenario 
given the additional information, i.e. the fact that the emission has 
such a large FWHM as compared to the b values of the absorption 
profiles and the near unity covering factor of the absorption gas. 
Therefore, the original conclusions of BOO still hold in which the 
absorption gas resides in a shell-like distribution, but we are now in 
a situation where the metallicity and the density of the absorption 
shells varies from shell-to-shell in the same source and also from 
source-to-source. 

U = -. — S , where Q = f kr- di/,ris the distance to the cloud from 

the ionising source and n is the particle density. 




0.001 0.01 0.1 1 



U 

Figure 6. The column ratios N civ / Nhi in photoionisation models as 
a function of the ionisation parameter U (represented by open diamonds 
joined by a solid line). The lower and upper curves correspond to gas with 
solar and 10 times solar metallicity, respectively. The dashed-line corre- 
sponds to models whose temperature has been increased in steps starting 
from the equilibrium value of 1 3000 K (open squares for T < 82000 K, 
filled triangles at higher T). The Nciv / Nhi column ratio actually de- 
creases when the temperature exceeds 10^ K. All models correspond to 
slabs for which the Nhi column is lO'^* '^^ cm"^ as in 0200+015 (ab- 
sorber #2 in Table 2). The dotted horizontal line represents the observed 
value of Nciv /Nhi ratio in 0200+015. 

6 DISCUSSION 

6.1 Origin of the absorbing shells 

In this section we speculate on the origin of the large absorption 
shells around HzRGs. 

6.1.1 Shells propagated by the passage of radio jets 

The evidence presented in this paper suggests that the absorbing 
shells surround the whole of the radio source and probably have a 
roughly spherical distribution. If this is the case then it is extremely 
unlikely that the shells are propagated by the passage of the radio 
emission for which we would expect the shells to be disrupted at 
radii less than the extent of the radio emission. The fact that we see 
the absorption systems in the majority of the smaller radio galaxies 
(projected linear size D < 50 kpc) but not in the larger sources 
(e.g. van Ojik et al. 1997) makes this scenario very unlikely. We 
also observe the absorption systems along the line-of-sight to the 
nucleus. Given that radio galaxies are orientated along the plane of 
the sky according to orientation based unified schemes, then this 
also rules out the scenario of a radio jet propagated shell. 

6.1.2 Starburst driven superwinds 

One explanation is that the shells originate at the centre of the 
galaxy and are ejected by superwinds, like those seen in Lyman 
Break galaxies (e.g. Steidel 2001). These winds will propagate 



8 Jarvis et al. 



roughly spherically, provide an explanation of the enriched nature 
of the shells and have no problem in explaining the number of 
shells and the diverse metallicities if the starburst driven winds are 
episodic. 

Evidence is mounting that powerful radio galaxies are formed 
at an early epoch and the triggering of the AGN activity is synchro- 
nized with a major episode of star formation (e.g. Archibald et al. 
2001; Willott et al. 2002; Baker et al. 2002). Therefore, it seems 
likely that starburst driven superwinds may be an important mech- 
anism for ejecting gas out of the central regions and the high dust 
concentrations inferred from submillimetre observations imply that 
starburst may easily be obscured at UV and optical wavelengths. 

Martin (1999) and Heckman et al. (2000) estimated the out- 
flow speeds of superwinds for galaxies of varying mass. They show 
that the outflow speed is largely independent of galaxy mass, imply- 
ing that the outflows preferentially escape less massive halos. High- 
redshift radio galaxies seem only to be found in the most massive 
elliptical galaxies at any redshift (e.g. Jarvis et al. 2001b), hence 
presumably the deepest potential wells. If this is the case then un- 
like the case for L* and sub-I/* galaxies the winds may not possess 
the required velocity to fully escape the deep potential well and 
subsequently enrich the intergalactic medium beyond the confines 
of the HzRG dark halo. Instead this enriched material resides at 
large radii, away from the central cusp, but still within the potential 
well of the galaxy. 

In the event of episodic superwinds, these shells would even- 
tually collide at large radii, which could both deplete or increase 
the metallicity of the most distant shell depending on the composi- 
tion of the wind. If this is the case then these shells will eventually 
fall back towards the centre of the shell in the absence of a repul- 
sive force. Indeed, with reference to Tables |l]and|, there is a slight 
tendency for the Lya absorbers in 0200-1-015 to have lower neu- 
tral column densities and b-parameters as the absorption redshift 
(Zshs) increases. If the difference between ^abs and the redshift of 
the emission line peak reflects outflow of the shell relative to the ra- 
dio galaxy, this trend could reflect evolution in the ionisation state 
and internal structure of the shells as they are ejected from the radio 
galaxy. 

This is in line with an initial major burst of star formation 
in which much of the gas is expelled, subsequently followed by 
smaller bursts in which the lower velocity, lower-density shells are 
expelled. At some stage in the evolution these shells would fall 
back under gravity and the shells would merge, but this process is 
probably very rare because the radio jet has already disrupted the 
outermost shells which are now seen in emission. However, it is 
difficult to reconcile the low metallicity shells in 0943-242 with a 
starburst scenario, as we would expect much higher enrichment if 
the shells originate in a starburst driven superwind. Although this 
is not the case for the high metallicity shell in 0200-1-015. 



6.1.3 Relics from a gas-rich merger 

Metal enriched shells of low-density gas could be a relic from a 
gas-rich merger in which there was an episode of star-formation, 
presumably induced by the merger. This scenario has been used to 
explain the vast HI shells in low-redshift elliptical galaxies (e.g. 
Sahu et al. 1996; Morganti et al. 1997). 

A number of numerical investigations predict that the outer 
regions of a galactic halo should be less metal rich than the central 
regions (e.g. Bekki 1998). This is in line with the case of 0943- 
242 but in opposition to the case of 0200-1-015 in which the main 



absorber is highly enriched and is comparable in metallicity to that 
expected in the narrow-emission line region towards the centre. 

An argument against the merger relic scenario is that the effi- 
ciency of mergers depends critically on the relative velocities, with 
low velocity mergers being the most efficient. This means that these 
efficient mergers would produce very low velocity dispersions in 
the absorbing gas. We might also expect the low density halo re- 
sulting from a gas-rich merger to have a more filamentary structure 
with a covering factor less than unity, which although reconcilable 
with our observations of a unity covering factor along our line-of- 
sight, is not the most plausible explanation. Therefore although a 
gas-rich merger remains a candidate for producing the extended 
low-metallicity haloes it has difficulties in producing the enriched 
shell we see in 0200-1-015 and probably needs a very specific merg- 
ing scenario. 



6.1.4 Gas-rich clusters and cooling flows 

It is now well established that many powerful radio galaxies reside 
in significant galaxies overdensities at all cosmic epochs (e.g. Hill 
& Lilly 1991; McLure & Dunlop 2001; Best 2001; Pentericci et 
al. 2000). Indeed, the discovery of ~ 20 Lya emitters around a 
z — 4.1 radio galaxy (Venemans et al. 2002) suggests that power- 
ful radio galaxies are ideal beacons for tracing proto-clusters in the 
early Universe. If this is ubiquitous for all powerful radio galax- 
ies then we might expect large reservoirs of gas around the radio 
galaxy. 

Within a cluster, cooling flows provide a mechanism for de- 
positing substantial amounts of cool gas and dust around the cen- 
tral galaxy (e.g. Edge 2001). The observational consequences of a 
radio source being triggered at the centre of such a cooling flow 
cluster at high redshift were discussed by Bremer et al. (1997). 
One of their clearest predictions concerns the presence of HI ab- 
sorption in the Lya emission in systems seen as radio galaxies (as 
opposed to radio-loud quasars); to see absorption from an ensem- 
ble of clouds requires a covering fraction of at least 50 per cent; 
little absorption is expected from clouds at radii of 10-100 kpc, as 
their covering fraction is low, but those between 10-20 kpc pro- 
vide a covering fraction of order unity and absorbing HI columns 
of 10^* - 10^1 cm ^ , as the surrounding gas away from the ion- 
ising beam should be similar to that of the interstellar medium be- 
fore the nucleus became active.. This would account for the ab- 
sorption in 0943-242 where the radio source has a projected linear 
size ~ 26 kpc; in 0200-I-0I5, however, the radio source is larger 
(~ 43 kpc) and has already engulfed these high-covering factor 
clouds, so the absorbers in this object may have been disrupted. 

However, the classical view of a cooling flow, from low- 
redshift studies is almost certainly not applicable at these high red- 
shifts. The main reason for this is that the environmental gas may 
not have had time to both settle towards hydrostatic equilibrium and 
also to be heated up to X-ray temperatures. In the low-redshift Uni- 
verse cooling flows also have an inclination to produce filaments 
when the gas cools (e.g. Fabian et al. 1981), which would result 
in a much reduced covering factor (^ 1; see e.g. Conselice, Gal- 
lagher & Wyse 2001). Therefore, a possibly more accurate view is 
one in which the pristine gas is collapsing and cooling in a more 
globally ordered fashion, i.e. as a spherical collapse. In this case 
the resulting gaseous distribution would be different from that in 
a low-redshift cooling flow and have the observed high covering 
factors. 
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6.2 A scenario: deposition of pristine gas and triggered star 
formation 

An age based scenario, which assumes that both the AGN activity 
and a major episode of star formation are synchronized, probably 
due to the collapse and cooling of (proto-)cluster gas (or a gas- 
rich merger), may explain the differences in the absorption shells of 
0200+015 and 0943-242. The projected linear size D, of the radio 
emission may be used as an age indicator for the radio source (e.g. 
Blundell, RawUngs & Willott 1999; Baker et al. 2002). In this case 
0200+015 (D ~ 43 kpc) is older than 0943-242 {D ~ 26 kpc). 
If this is the case then there is much more time for the absorption 
shells to have evolved in 0200+015. An indication of this comes 
from the interaction of the radio emitting plasma in the extended 
component in 0200+015, which looks to be disrupting the absorp- 
tion shells and subsequently producing additional emission. 

If starburst or quasar driven winds precede the passage of the 
radio activity then a pristine halo at large radii will be enriched 
by the winds before the radio plasma propagates through the ex- 
tended regions. Indeed, in an investigation of the metallicity of the 
broad-line region in quasars, Hamann et al. (2002) compared the 
UV emission line ratios in a sample of radio-quiet quasars with 
theoretical models. They found the broad-line region gas was most 
likely to be of high metallicity {Z > Zq) and was a consequence 
of rapid star formation triggered at or before the time of the quasar 
activity. Ejection of the gas via a quasar or starburst driven outflow 
would obviously enrich the surrounding halo and thus offer an ex- 
planation of the highly enriched nature of absorber #2 in 0200+015. 

In the case of 0943-242 any starburst driven wind associated 
with the activation of the nucleus may not have had time to enrich 
the outer halo. Thus the explanation of Binette et al. (2000) still 
holds for this source, in which the shell has a very low-metallicity. 

This would also explain the variety of metalUcities observed 
in these shells, as some of the shells will have been enriched, while 
the winds have not yet reached the others. This naturally fits in with 
the absorption lines in 0200+015, in which the the blueshifted (ab- 
sorber #2) shell has a high metallicity, whereas absorbers #3 and 
#4 he redward of the Lya emission line and are probably infalling 
components from the (proto-)cluster environment, in which there 
has been little or no star formation to enrich the shells. This is also 
in Une with the work of Overzier et al. (2001) in which the metallic- 
ity gradient in the emission-line gas in the gaseous halo surrounding 
MRC2104-242 may also arise from the same mechanism. 

The fate of the shells is uncertain, it is not implausible that 
shredding by the passage of the radio jet (e.g. Bremer et al. 1997) 
will work to condense the shells into distinct clouds that may re- 
semble the knots of emission-line gas observed in many larger ra- 
dio sources. This naturally fits in with correlation between the ex- 
tent of the emission-line gas and the extent of the radio emission 
(e.g. Jarvis et al. 2001a), with the radio emitting plasma propagat- 
ing through the low-density gas which consequently increases the 
density of this gas resulting in detectable line emission. 



6.3 The link to quasars 

Using a complete sample of radio-loud quasars from the Molon- 
glo Quasar Catalogue (Kapahi et al. 1998), Baker et al. (2002) used 
both the Hubble Space Telescope with STIS and ground based tele- 
scopes to investigate the associated C IV absorption in two redshift 
intervals. They found that the absorbing systems in these objects 
have similar velocities widths and offsets to the two radio galaxies 
investigated in this paper, the gas is mixed with a large amount of 



dust and that the gas and dust dissipate over the lifetime of the ra- 
dio source. Thus, the main conclusions from this study fit in easily 
with our results, in line with Unified schemes. 

The advantage that radio galaxies hold over quasars is that, ac- 
cording to unified schemes, the radio galaxies are orientated along 
the plane of the sky and so the deprojection to the true linear size of 
the radio source is trivial. Therefore, age effects may become more 
apparent in a large sample of radio galaxies with higher signal- 
to-noise observations. This should easily be achievable with cur- 
rent instrumentation as we have shown that spectra of resolution 
of ~ 2A at these redshifts is sufficient to probe the absorption 
shells to the required detail allowing the integration time to be sig- 
nificantly reduced. Combined with a similar sample of radio-loud 
quasars, in which the reddening properties may be explored, we 
should be able to track the evolution of both the absorption- and 
emission-line gas and the intrinsic dust properties in all radio-loud 
AGN. 



6.4 Further observations 

Obviously further observations of the sources studied in the paper 
will enhance our understanding of the origin and fate of the gas sur- 
rounding HzRGs. Indeed, deep spectroscopy around other resonant 
lines may provide critical information regarding the state of the ab- 
sorbing shells. The most obvious emission line to use for this study 
would be the Mgll AA2795,2802 doublet. The absorption equiva- 
lent width ratio of C iv/Mgll is frequently used to classify quasar 
absorbers into high or low excitation (e.g. Bergeron et al. 1994) 
and will also provide further information regarding the metalUcity 
of the gas. 

0200+015 is also at a redshift in which the Ha emission line 
is redshifted into the near-infrared if— band ttansmission window. 
As Ha is not a resonant line the true velocity structure of the Lya 
line should be replicated in Ha, thus allowing a detailed study of 
the actual kinematics of the emission-line gas. The Lya/Ha ratio 
is also sensitive to reddening by dust, and if the Lya photons es- 
cape from the clouds via resonant scattering from HI, the Lya line 
should have significantly broader wings than that of Ha. 

Finally, spectroscopic observations around various resonant 
line of a large sample of high-redshift radio galaxies spanning pro- 
jected Unear sizes from D ~l-2 kpc up to the classical doubles 
with D > 100 kpc are warranted. These will shed further light on 
the distribution, density and metallicity of gas in the environments 
of the most massive galaxies at high redshift, and provide clues as 
to the origin and fate of these large gaseous halos, and may conse- 
quently provide important clues to galaxy formation in general and 
the important role that feedback may play. 



7 CONCLUSIONS 

We have obtained high-resolution spectra of two high-redshift ra- 
dio galaxies to determine the structure of the absorption gas which 
appears to be ubiquitous in small (D < 50 kpc) radio sources. The 
main conclusions of this study can be summarized as follows: 

• The higher resolution observations do not uncover significant 
structure in the absorbers of both Lya and Crv in 0943-242, how- 
ever the A^Hi ~ 10^® cm"^ absorber in 0200+015 split into two 
absorber of A^hi ^ 4 x 10^'* cm~^. The absorption gas present 
in our two chosen high-redshift radio sources is most Ukely dis- 
tributed in shell-like structures which encompass the radio source. 
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• The metallicity of these absorption shells can vary between 
sources and also between the shells in the same sources, implying 
that the shells are enriched by a secondary process. 

• We speculate that the most plausible origin of the absorbing 
shells is two fold. The pristine shells are the relics from a gas-rich 
merger or (proto-)cluster gas which has cooled and collapsed to- 
wards the centre of the dark matter halo. Whereas the enrichment 
comes at a later stage when the AGN is triggered simultaneously 
with a major episode of star formation, presumably due to a large 
inflow of gas into the nucleus. This major burst of star formation 
initiates starburst driven superwinds which drive the gas toward the 
outer halo of the galaxy and may explain the very high metallicity 
of the absorption shells in 0200-1-015. 

• The passage of the radio emitting plasma disrupts and frag- 
ments the absorption shells in the larger, and thus older radio 
sources. Evidence of this is seen in 0200-1-015, where the absorption 
screen along our line-of-sight to the radio lobe appears to be more 
fragmented than the absorbing screen between us and the nuclear 
emission. 

• We have found direct evidence that the narrow-emission lines 
in high-redshift radio sources may be outflowing from the central 
regions. This is in agreement with previous studies of the narrow- 
emission line dynamics in nearby Seyfert galaxies and is in line 
with orientation based Unified Schemes in which the outflow from 
the far side of the source is obscured by dust. 
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